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Abstract 

Based on the kinetic theory, the non-Abelian medium property of hot 
Quark-Gluon Plasma is investigated. The nonlinearity of the plasma comes 
from two aspects: The nonlinear wave-wave interaction and self-interaction 
of color field. The non-Abelian color permittivity is obtained by expanding 
the kinetic equations to third order. As an application, the nonlinear Landau 
damping rate and the nonlinear eigenfrequency shift are calculated in the 
longwave length limit. 
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I. INTRODUCTION 

Looking for Quark-Gluon Plasma (QGP) as an important aim of relativistic heavy ion 
collisions has inspired wide interest and stimulated extensive studies 0^. An urgent and 
important subject is to investigate medium effects in the hot plasma. The high transverse 
momentum particles and jets created by the hard scattering processes, which propagate 
through the hot and dense matter produced by the soft processes, can be used as a clear and 
observable probe of QGP 00 . Theoretically, to study the behavior of high energy partons 
through QGP one can construct a kind of response theory for QGP to external current, 
which must be non-Abelian and hence nonlinear 0. The non-Abelian permittivity, as a 
characteristic quantity of QGP medium, must be known. 

The kinetic theory is a powerful tool to investigate the medium effects. The QGP kinetic 
theory which is fundamentally capable of dealing with the non-equilibrium phenomena has 
been formulated (see, for example. Ref. |]^), and the Abelian-like permittivity of QGP in 
the linear approximation has been given out 0. However, in the linear approximation, the 
property of QGP is quite similar to QED plasma. 

The great difference between the QED and QGP plasma is that the former is an Abelian 
medium while the latter is a non-Abelian one. It is not sufficient to take only the Abelian 
dominance approximation or the linear approximation in solving QGP kinetic equations. In 
order to reflect the essential characteristic of QGP, one must go beyond the linear approx¬ 
imation and take into account the nonlinear wave-wave interaction and the self-interaction 
of color held, as pointed out by Ref. |TB| . 

In this paper, we aim at discussing the nonlinear medium effects in QGP based on the 
kinetic theory of chromodynamics. By using the collisionless Vlasov kinetic equations, the 
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induced color currents up to third order are obtained. Then from the color held equation 
including both the nonlinear wave-wave interaction and the self-coupling term, the color 
permittivity of QGP is calculated. As an application, the nonlinear eigenfrequency shift 
and the nonlinear damping rate are calculated for the mode k = 0. 

The paper is organized in the following. In Section we review the kinetic equations of 
QGP and discuss the expansion of the distribution function in the color mean held. Section 
H is devoted to obtaining the hrst order, second order and third order color currents by 
solving the kinetic equations. In Section |V|, the general expression of nonlinear permittivity 
is obtained by substituting the obtained color currents into the color held equation. The 
nonlinear Landau damping rate and the nonlinear eigenfrequency shift are obtained from 
the color permittivity in Section 0. 

The natural units and the Minkovski metric are used through the paper, i.e., h = c = 
ks = l and diag{g^,^) = diag{g^^'') = (1, -1, -1, -1). 


II. KINETIC EQUATIONS FOR FLUCTUATIONS 
A. Mean Field and Kinetic Eqnations 

The color mean held equation in QGP is 

D,F^’^{x)=f{x), ( 2 . 1 ) 


with the color current 

fix) = I / (^2TiyE “ Q-ix,p)) + ‘2itafabcGbcix,p)] (2.2) 

and the non-Abelian mean held stress tensor 

Fff) = d^Afx) - a^A“(x) + gfabcAlix)Al{x). (2.3) 

Here, t°‘,a = 1, 2, • • ■, N^ — 1, are the SUiNc) group generators in fundamental representation 
and fabc are the corresponding structure constants. 

The quark and anti-quark distribution functions Q±{x,p) are Nc x matrices in 
color space, and behave as Qix,p) —>• Uix,p)Qix,p)Wix,p) under the transformation 
operator Uix,p). The Gluon distribution function Gix,p) transforms as Gix,p) 
M{x,p)Gix,p)M'^{x,p) with Mab = Tr[t“U(x,p)t^t/'^(a;,p)]. 

In the collisionless semi-classical limit with spin being neglected, the distribution func¬ 
tions Q±{x,p) and G{x,p) satisfy the following transport equations §’i 

p^Df,Q±{x,p) ± ^p^^{Ff,„{x),Qfx,p)} = 0 , 

2 dpu 

p>^D^Gix,p) + |p^—{JQ,^(a;),G(a;,p)} = 0, (2.4) 

where and are the covariant derivatives dehned as 
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Df^ dfM ig[A^, ' ‘ 

D^ = d^-ig[A^,---]. (2.5) 

The color mean field A^{x) and Afj, and the field stress tensors F^i, and have the following 
Lie-Algebra relations 


A^{x) = AC(a;)r, F^%x) = Fr{x)F, 

Aab{x) = -ifabcA^ix), Kb{x) = -ifabcF^''{x). ( 2 . 6 ) 


The Eqs. (p.l|) and (p^ ) form a closed set of coupled equations. It is easy to see that the 
nonlinearity in QGP results from two aspects. One is associated with the nonlinear kinetic 
equations and the other one is related to the self-interaction term on the left hand 

side of mean field equation (p.l|) . 

An usually used approximation to solve the above coupled equations is the linearization. 
Linearizing the left hand side of mean field equation (|2.4|) , i.e., neglecting the self-interaction 
term, was considered in Refs. [0,^], and linearizing the kinetic equations was discussed by 
Blaizot et ah in Ref. . What we will do in this paper is to discuss the nonlinear effects 
in both the field and kinetic equations. 


B. Equations of Fluctuations 


For convenience to analyze the wave-wave interaction of fluctuations in QGP, the mean 
fields and the distribution functions can be divided into regular parts and fluctuation 
parts ili’Ei0iil. The quark, anti-quark and gluon distribution functions are decomposed 
as 


Q±{x,p) = Q±(x,p) -f g±(a;,p), 

G{x,p) =G{x,p)F g{x,p), (2.7) 

and in the meantime, the field A and current are decomposed as 

= A“(x) + Af{x), j^{x) = j“(x) + j^^{x). (2.8) 

Here the quantities having a bar represent the ensemble average of the corresponding quan¬ 
tities, and by definition the mean values of the fluctuations vanish 

{q±ix)) = 0, {g{x)) = 0, 

(A'^ix)) = 0, {fix)) = 0. (2.9) 


With the above decomposition, one can obtain the kinetic equations for the fluctuations 
of distribution functions by subtracting the kinetic equations for the regular distributions 
from the average of the kinetic equations (|2.4|) 0 


p^df,q±ix,p) = ig{[Al{x),q±{x,p)] - ([Aj(x), g±(x,p)])) ^ -gp^iF^^if), 


dQ±ix,p) 

dpu 


} 






dpv 
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±^^vY{lAl(x),Al(x)] - {[A^^(x),Al(x)]},?^^^}, 
lfd^g(x,p) = ig[lAl[x),g(x,p)] - (\Al(x),g(x,p)])) - \gp>‘{Tl^x(x), } 


dpu 




}-({^J„.M,%^})) 


■ii9V({[9lJ(a:),9lJ'(x)],^^} - {{IAHx),AHx)],?^^}}) 
ii9V{K(9^).9lJW] - ([9lJW,9ir(x)|), 


( 2 . 10 ) 


where the index L denote the linear term of in A^. 

The fluctuation part of the color held satishes the following Yang-Mills equation 

a.Fp" + = -ill + 21 ^"! - a^"])) 

+»9([Y. Y'I - (lY. Y'D) + al\Al, [A^^A^-]] - {[Al, ^l^-]])), (2,11) 

where the current huctuation is also divided into a linear term and a nonlinear term 
jjA- By expanding the huctuation of the distribution functions in terms of the huctuations 
of the color held AA 


00 

q±{x,p) = Y, q±\x,p), 

n=l 

00 

9{x,p) = 

n=l 

the nonlinear color current can be expressed as 

CXD 

71=2 


with 


d^P ^ 

{2n)^p^^ 


tr (f“ {q+ \x,p) - + Tr{T^g^^\x,p))] , 

n = 2, 3,4, • • ■. 


( 2 . 12 ) 


(2.13) 


(2.14) 


Substituting the above expansion of huctuations in Aj^ into the kinetic equations ( p.lOD , 

one obtains the hierarchies of kinetic equations for the distribution functions and 

By using the regular distributions and regular held as input, which can be calculated 
from the corresponding kinetic equations and mean held equatiom one can in principle solve 
the above kinetic and mean held equations for the huctuations oBMEJ. For simplihcation, 
we will consider in the following the huctuations above the colorless equilibrium state. We 
take y4 = 0 and Fermi-Dirac or Bose-Einstein distribution function for Q± and G 
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g/3p-M _ X ’ 


(2.15) 


Q± 


2iV/ ^ 

el3p-u ^ 1 ’ 


where u = (1, 0, 0, 0) is the local four-velocity of the plasma, Nf is the the number of flavors 
for quarks. Since we are interested in the fluctuations only, the index T for the held and 
current will be omitted from now on. 

It’s convenient to discuss the huctuations in momentum space through the Fourier trans¬ 
formation. 


r 

q±{x,p) = J (2.16) 

We choose temporary axis gauge, i.e., Aq = 0, and consider only the longitudinal held 
excitation, A’‘{k) = with iF = |k|. 


III. LINEAR AND NONLINEAR COLOR CURRENTS 
A. Linear Approximation 

In order to investigate the nonlinear property of QGP, one should solve the kinetic 
hierarchies for the distributions q±'^ and g^') order by order. Let’s hrst review the linear 
approximation which has been discussed in many previous works. The dispersion relation 
obtained here will be used to perform the relevant integrals in the end of this paper. 

gFrom the kinetic equations for the hrst order huctuations of the distribution functions, 
it is easy to write down 


q±\k,p)= T 


k^p^A^k) 

p ■ k + ip^e~^ 


gW{k,p) = 


k^p'^A^{k) 
p ■ k + ip^e~^ 


dQ±{p°) 

dpO 

dG{p^) 

dpO 


(3.1) 


Substituting these hrst order huctuations into (p^.l4| ), one obtains the hrst order color current 
by noting TrT“T^ = Nc6°‘^. Then from the mean held equation to the hrst order, one has 
the nontrivial solution 


— k)A{k) = 0, 

with the well-known linear dispersion relation 

eL{uJ,k) = 1 + ^[1 - ^(Inl ^^^ l -i7Te{K-u))] = 0, 

where 


LJp — 


'(IV/+ 21V.) 


18 


9^2 


(3.2) 


(3.3) 


(3.4) 
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is the plasma frequency. In the longwave length region, the dispersion relation is reduced to 

o;^ = (Up + (3.5) 

It is necessary to note that the linear dispersion relation (|3.3|) means that the eigenwaves 
in QGP are always timelike, they can not exchange energy directly with the particles in 
the plasma. Therefore in the linear approximation, the eigenwaves in QGP do not undergo 
damping fl. 

To obtain the nonlinear effects in the medium, one should solve the kinetic equations to 
nonlinear orders, i.e., to go beyond the linear order approximation. 


B. Second Order and Third Order Color Cnrrents 


In this subsection, the second order and third order color currents will be calculated. 
The following shorthands 


^ = y 6{k — ki — k2)dkidk2- 


k=ki+k2 


(2vr)^’ 

= J 6{k — ki — k2 — ks)dkidk2dks 


k=ki-\-k2-\-k3 


(2nf' 


(3.6) 

(3.7) 


are used for notational simplihcation. 

The second order kinetic equation for the quark distribution is the following 

- W^k^q^+\k,p) = igp^ {[Mh),q+\k2,p)] - {[Ai{ki,p),q^+\k2,p 

ki+k2=k 


+^9if E 


ki+k2=k 


dpi 


dq+\k2,p) 

dpi 


}) 


^ k=ki+k2 \ OPi, 

+>‘ E ({[.4.(fci)..4,(fc,)]. ^yy;^'^d - ({|.4.(fci).^i(*^2)i. 

/c=/ci+fc2+^3 ^ ' 

The kinetic equations for the second order fluctuations of the anti-quark and gluon dis¬ 
tribution functions can be written out similarly to ( p.8|) with opposite signs for the terms 
related to {■ ■ ■} for anti-quarks and with the replacement of q, A and by G, A and 
for gluons. 

Inserting the second order distribution functions into Eq. (|2.14|) , one obtains the second 
order color current 


J 


T{2)ka _ 


d^p p^pip^ 




^9 fabc J - k + ip^’e p ■ kq + ip^s 


- {Atik^AjihA) . (3.9) 


(3.8) 
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(3.10) 


with the effective equilibrium distribution function 

The third order color current can be obtained similarly 

jT{3)la^k^ = 

E {At(.h)A'‘(M)Al(k2} - AKk.XAHkOAUh)} - {AtikMlikMUh)}] . (3.11) 

k=ki+k2+k3 

with 


Y^abdeijkl _ rabc fcde-sp{I)ijkl , ^ab^de-sp(II)ijkl , iabcjcde-sp{III)ijkl 

^kkik2k3 — J J ^kkik2k3 "V 0 0 ^kkik2k3 “ ^kkik2k3 
^gatgde ^ g^gte _g_ S'^S'^)SilXf‘, 

d?p p^pip^p^ 1 u!2dp„'D, 


^{I)ijkl _ 4 

^kkik2k3 y 


■po^eq 


y{II)ijkl ^ 

^kk\k2k3 jy ‘^kkik2k3 

9^ 


(27r)3p° p ■ k + ip^e p ■ (fci + ^2) + ip^e p ■ k2 + ip^e ’ 

{III)ijkl 


d^p p^x^^{k3,p) d 


X^’'{ki,p) d 


UJ2dpeq 


2Nc J (27r)3p° p ■ k + ip^e dpf^ \p ■ {ki + ^ 2 ) + dp^ p ■ k2 + ip^e' 


y{VI)ijkl 

and 


r 
’ 2 


d^p p^x*^(fc3,p) d 


X^’'{ki,p) d 




(27r)^p° p ■ k + ip^e dp>^ \p ■ {ki + ^ 2 ) + ip^e dp’-' p ■ k2 + ip^e' 


(3.12) 


X^’'{k,p) = {p-k)g^^^-p^k’', (3.13) 

where dbce are the SU{Nc) symmetric structure constants. From Eqs. (^^) and (|3.11|) , one 
can see that the second and third order color currents come from the response of particles 
to the secondary waves resulting from the nonlinear interaction of eigenwaves. 

Before making further calculations, one should note that on the right hand side of ( p.ll ), 
the term with purely non-Abelian, while the others are Abelian terms. Since the 

soft excitation momentum k ~ gT and the particle momentum p ^ T i’i, one can easily 
conclude that the leading non-Abelian term is proportional to while the Abelian terms 
are proportional to the higher orders of the coupling constant g: ~ ^kkikSf ~ 9^ 

and g^ and ~ g\ 


IV. NONLINEAR PERMITTIVITY 

Now we turn to calculate the non-Abelian permittivity from the mean held equation to 
the third order 

- kpk^^A\k)+kjk^A^{k) + = 

,/ E (l3l,(fci).lV(fa),4-(A.3)|]-(14,(fa),|4J(fe),4‘(A.3)]]))-iti(fc), (4.1) 

k=ki+k2+k3 
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with 




(4.2) 


/^From ([4. 11 ) , one can see that the left hand side is Abelian-like, while the right hand side 
is non-Abelian which results from two different origins. The hrst term on the right hand 
side of (^4.1|) with the summation represents the color self-interaction. Some of our previous 
works [l^J^ neglected this nonlinearity. The second non-Abelian property in (|4.1|) is from 
the nonlinear current jj*. This nonlinear wave-wave interaction is not taken into account 
in Refs. [0|,^,pT| . Here we will discuss the solution of the held equation by considering the 
color self-interaction and the nonlinear currents simultaneously. 

Substituting and into the held equation (^4.1|) , one gets the relation 




K 


3+a £abc 


igH-f 


E 

k=k\+k2 


(Pp (p ■ k) (p ■ ki) (p ■ ka) 


UJ2' 


dVeq 

dpo 


-gH<^ 


(27r)3p0 K Ki 
X {A\k,)A'^{k2) - {A\kM\k2))) 

Y- \ 2 f (p-ki)(p-k2)(p-k3)(p-k) 

k=ki+k2+k^ - 


K2 p ■ k + ip^0+ p ■ k2 + 
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(27r)^p° 


K^K2K^K 


p ■ k + ip^e 


X 




po^eq ^ (ki ■ ka) (k ■ ks) 


p- {ki + /ca) + p-k2 + ip^e ' K1K2 KK3 
X [A\k^)A\k{)A\^2) - A\k^){A\k{)A\k2)) - {A\ks)A<^{k,)A%k2))) • (4.3) 

Obviously the hrst and second integrals rehect the three-wave process {k ^ ki + fca) and the 
four-wave {k ^ ki + k 2 + k^) process, respectively. 

It is well known that it is difficult to solve the general Yang-Mills equation because of 
the color held self-coupling term. However, if the phase of excitations in QGP is random 
(Gauss-like distributed), the correlation of even random ouantities can be approximately 
expressed as the product of correlations of two quantities @§ 1 , and the solution of the held 
equation can be greatly simplihed. Concretely, after multiplying both sides of (^^) with 
A^{k'), taking averse with respect to the Gauss-like random phase, and then using the 
following relations £3 

= {A^A^)^,^, {A\k)A^{k')) = (27r)^/f 5(fc + fc'), (4.4) 

{A^ik^)A\k 2 )A^{ks)A\k,)) = (27r)8h(fci + ^2 + ^3 + k,) x [6{k, + k 2 )lfjg 

+6{k, + + 5{k, + , (4.5) 


the equation ( |4.3| ) is reduced to 


+ £*ad) dk^ — 0 , 


(4.6) 


where is the linear permittivity given in (|3.3|) , and the nonlinear permittivity contains 
two parts 
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with 


£ 


* 

ad 



+ S 


5”+ 
ad 5 


s* ^ 9 

^ad ^2 


cabc £ced jbe 


R* ^ ^ 
^ad 2 


rabc rc 

J J -‘fci 

'■ (Ppdp^Ve 


( ^■^1 

[kk. 


abc 4 ^cde jbe i j^adc j^cbejbe\ 


{nf 


fci 


d/ci 


-’po^eq rabc 


(27r)^p° 

1 


r 


'pki\ /p-k' 






p ■ /c + ip^e 


UJ 


X 


p ■ (fc — fci) + ip°e \p ■ k + ip^e 


fcedjbe | fcdejbe\ dk\ ^ ^ 

J -'fci T 'j ; ^ITJ -'fci /r, \A- 

^ p ■ ki + ip^e ^ J {271)* 


Here and correspond to the first and second term on the right hand side of ( |4.1|) , 
respectively. They reflect the purely non-Abelian medium effects of QGP. From Eq. o, 
the nonlinear permittivity is a matrix in color space because of the non-Abelian SU{Nc) 
structure constants. Its value depends on the correlation strength = {A^A^)k and the 
dispersion relation a;(k) of the collective excitation. 

For simplihcation, we consider only the correlation between the same colors and take a 
diagonal matrix in the color space 


— 2 ^d{u) — CUk) + d{u) + U7]i))Ik.Sbc-, 


TT 


U 


(4.8) 


where Jk is the intensity of th^TOirelation with frequency Wk. For fluctuations at thermal 
level, one can take Jk = dvrT 
color permittivity are further reduced to 


Then the nonlinear contributions efj and to the 




2g^N, 

a;2(27r) 


,k ■ ki 

KKi 


^d^ku 


o;, 


ki 




2g'^Nc f (PpdpVeq /p-ki\^ /P'k' 




(27r)^p° 

1 




K 


X 


UJq 


p ■ k + ip^e 
1 


^ki 


,3 


^p ■ {k — ki) + ip^ep ■ k + ip^e p ■ {k + ki) + ip^ep ■ ki + ip^e J a; 


(4.9) 


ki 


p.(k-ki)+ipO€ expression of e^* have dif- 


The four factors p.l^_^_^p 0 ^y p.k-i^^ip0^jp.(^k+ki)+ip°€ ^ - 

ferent contributions to the integral. The hrst three factors have only real contributions 
because the eigenwaves are always timelike, while the last one has both real and imaginary 
contributions to due to the well-know formula 


= P- 


TT 


p ■ {k — ki) + ip^e p ■ {k — ki) 


- i-r(5(a;k - Wki - V ■ (k - ki)), 

pU 


(4.10) 


where v = p/p°. The imaginary part of e^* indicates that the secondary waves resulting from 
the nonlinear interactions can exchange energy with the plasma particles, which corresponds 
to the nonlinear Landau damping. 

To perform the above integrals over ki and p, we can choose spherical polar coordinates 
and take the polar axis of ki in the direction of k and the polar axis of p in the direction of 
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ki. Considering the soft excitation in QGP 0, the np and low limits of Ki can be chosen 
to be gT and g'^T, respectively. 

By using Eq. ( [4 .101) and the integral 


{Nf + 2N^)7 i^T^ 


PodpfyD^qdpQ 


(4.11) 


the nonlinear color permittivities and are finally written as one dimensional integrals 






3(27r)4a;2 Jg^r 0 ;^^ 


2 KldKi, 


Re{e^*) = 


8g^N^7r^T^{Nf + 2N^ 


9(27r)"‘a;^ 
r9T KfdKi 2 
Ig^T CUki Kl 


(qT KfdKi 2(a;ki “ 


Ig'^T 10- 


1 + 


a;. 


ki 


‘2uj-pK\ 


In 


Kl 

(Uki — Ki 


1 - 


Ct^ki 


l^ki + Ki 


2Ki 

(^ki ~ Wp) 


In 


^ki 


— l^ki + Kl 


2u}pK\ 


^ki ^p 


^^ki ~ + Kl 




NcTT^T\Nf + 2Nc) [9T 1 


—dKi{uJk^ 

g^T Kl 


Upf — . 


9(27r)4a;^ Jg^T Ki ' ‘ ‘ 

The contributions of e^* and e^* can be seen quantitatively by making a rough numerical 
estimate. Using the dispersion relation (|3.5|) in the longwave length limit, the total nonlinear 
permittivity is 


(4.12) 


+ Re{e^*) + zlm{e^*) ^ -0.98c/ - 1.91^ + 0.61^^ = -2.89g + O.Ol^z. (4.13) 


V. CONCLUSION 


Based on the chromodynamic kinetic and mean held equations, the nonlinearity in QGP 
has been analyzed. One nonlinearity comes from the response of particles to the secondary 
waves resulting from the nonlinear interaction, and the other nonlinearity is due to the non- 
Abelian self-interaction in mean held equation. The non-Abelian permittivity due to the two 
nonlinearities has been calculated by considering the second and third order color currents. 
/^From its general expression, it depends on the dispersion modes uj\^ and the correlation 
strength of color held. The non-Abelian property is rehected in the presence of the SU{Nc) 
structure constants fate- With the given color permittivity, it is possible to further analyze 
the non-Abelian response of QGP to external currents 0. 

As an application of the nonlinear color permittivity ( 4.13|) , we can obtain the nonlinear 
Landau damping rate 0 

= ( 6 . 1 ) 


7 = 


deL I 

diO 


and the nonlinear eigenfrequency shift 


Ao; = 


Re{e* 


Pel I 

duj 


= O.SSjT, 


(5.2) 


which characterizes the nonlinear ehect in QGP too. 
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